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We study initial-state radiation events, e + e~ — > yi SR n + TT~ J/tp , with data collected with the BABAR 
detector. We observe an accumulation of events near 4.26 GeV/c 2 in the invariant-mass spectrum 
of ir + n~ J/tp . Fits to the mass spectrum indicate that a broad resonance with a mass of about 
4.26 GeV/c 2 is required to describe the observed structure. The presence of additional narrow reso- 
nances cannot be excluded. The fitted width of the broad resonance is 50 to 90MeV/c 2 , depending 
on the fit hypothesis. 



PACS numbers: 14.40.Gx 13. 66. Be, 13.25.Gv 

Recent observations of the A (3872), decaying into 
ir+ir-Jhb 0, Q, H 0, and the V(3940), decaying into 
u> J/tp 5| , have renewed experimental interest in char- 
monium spectroscopy. We have previously reported a 
search for direct A(3872) production in e + e~ annihi- 
lation through initial-state radiation (ISR): e + e _ — > 
7i SB A 6]. No signal is observed, suggesting that the 
A (3872) is not a 1"~ state, just as expected for a nar- 
row state well above the DD threshold. In this Letter 
we present a study of the e + e _ — > 7 IfiI! 7r + 7r~ J/tp process 
across the charmonium mass range. 

We use data collected with the BABAR detector at 
the SLAC PEP-II asymmetric-energy e + e~ storage ring. 
These data represent an integrated luminosity of 211 fb _1 
collected at ^/s = 10.58 GeV, near the peak of the T(4S) 
resonance, plus 22fb _1 collected approximately 40MeV 
below this energy. 

Charged-particle momenta are measured in a tracking 
system consisting of a five-layer double-sided silicon ver- 
tex tracker (SVT) and a 40-layer central drift chamber 
(DCH), both situated in a 1.5-T axial magnetic field. 
An internally reflecting ring-imaging Cherenkov detec- 
tor (DIRC) with quartz bar radiators provides charged- 
particle identification. A Csl electromagnetic calorime- 
ter (EMC) is used to detect and identify photons and 
electrons, while muons are identified in the instrumented 
magnetic flux return system (IFR). 

Electron candidates are identified by the ratio of the 
shower energy deposited in the EMC to the momentum, 
the shower shape, the specific ionization in the DCH, and 
the Cherenkov angle measured by the DIRC. Muons are 
identified by the depth of penetration into the IFR, the 
IFR cluster geometry, and the energy deposited in the 
EMC. Pion candidates are selected based on a likelihood 
calculated from the specific ionization in the DCH and 
SVT, and the Cherenkov angle measured in the DIRC. 
Photon candidates are identified with clusters in the 
EMC that have a shape consistent with an electromag- 
netic shower but without an associated charged track. 

A candidate J/tp meson is reconstructed via its decay 



to e + e~ or The lepton tracks must be well recon- 

structed, and at least one must be identified as an elec- 
tron or a muon. An algorithm to associate and combine 
the energy from bremsstrahlung photons with nearby 
electron tracks is used when forming a J/ip — > e + e~ can- 
didate. An e + e~ (fi + [j,~) pair with an invariant mass 
within tH (-11) MeV/c 2 of the nominal J/tp mass is taken 
as a J/tp candidate and is combined with a pair of oppo- 
sitely charged tracks that are identified as pions. 

Following an observation of an enhancement in the 
ir + ii~J/tp mass spectrum during an earlier search for 
ISR A" (3872) production in a 124fb _1 subsample of the 
available data, we chose to exclude the mass region 
from 4.2 to 4.4 GeV/ c 2 from consideration during opti- 
mization of the selection criteria with the full sample 
to avoid the introduction of statistical or other biases 
in the analysis of this region. Radiative production of 
the tp(2S) serves as a clean benchmark process @ for 
a data-driven optimization. Selection criteria are cho- 
sen to maximize A/ (3/2 + VB) 0], where N is the 
total number of 7 ISR V'(25'), tp(2S) — > ir + ir~J/tp candi- 
dates in the 20MeV/c 2 7r + 7r~ J/tp mass range that brack- 
ets the tp(2S) mass, and B is the number of events in 
the TT + TT~J/tp mass regions [3.8, 4.2] GeV/c 2 and [4.4, 
4.8] GeV/c 2 , scaled to the width of the originally observed 
peak. Simulated ISR events are validated with the tp(2S) 
data and are used to extrapolate the selection criteria to 
the excluded mass region as appropriate for small kine- 
matic differences due to the higher mass. 

Radiative e + e~ — > 7 KR 7r + 7r _ J/tp events are character- 
ized by a small mass recoiling against the tt + tt~ J/tp sys- 
tem and by low missing transverse momentum. These 
properties are reflected in (1), (2), and (3) of the se- 
lection criteria: (1) there must be no additional well- 
reconstructed charged tracks in the event; (2) the trans- 
verse component of the visible momentum in the e + e~ 
center-of-mass frame, including the ISR photon when 
it is reconstructed, must be less than 2.5GeV/c; (3) 
the inferred value of the square of the mass recoil- 
ing against the 7r + 7r~ J/tp combination (w|j ec ) must be 
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within [-1.02, +3.27] GeV 2 /c 4 for J/0 -> e+e~ candi- 
dates and [-1.06, +1.25] GeV 2 /c 4 for J/0 -> can- 
didates; (4) cos0£, where 6^ is the angle between the 
£ + momentum in the J/ip rest frame and the J/0 mo- 
mentum in the e + e~ center-of-mass frame, must satisfy 
| cos ^ | < 0.90. In addition, (5) for the e + e~ mode, 
cosOk, where 9^ is the angle between the n~ momen- 
tum and the J/0 momentum in the tt + tt~ rest frame, is 
required to be less than 0.90 to reject background from 
misidentified low momentum e~ in the forward region of 
the detector. We do not require the ISR photon to be 
detected in the EMC since it is produced preferentially 
along the beam direction. 

Candidate ir + ir~£ + £~ tracks are refitted, constrained 
to a common vertex, while the lepton pair is kine- 
matically constrained to the J/ip mass. The resulting 
7r + 7r~ J/0 mass-resolution function is well-described by a 
Cauchy distribution 0] with a full width at half max- 
imum of 4.2MeV/c 2 for the ip(2S) and 5.3MeV/c 2 at 
4.3GeV/c 2 . 

The it + tt~ J/ip invariant mass spectrum for candidates 
passing all criteria is shown in Fig.^ as points with error 
bars. Events that have an e + e~ (/i + /i _ ) mass in the J/ip 
sidebands [2.76, 2.95] or [3.18, 3.25] ([2.93, 3.01] or [3.18, 
3.25]) GeV/c 2 but pass all the other selection criteria are 
represented by the shaded histogram after being scaled 
by the ratio of the widths of the J/0 mass window and 
sideband regions. An enhancement near 4.26GeV/c 2 is 
clearly observed; no other structures are evident at the 
masses of the J PC = 1 charmonium states, i.e., the 
0(4040), 0(4160), and 0(4415) [n|, or the AT(3872). The 
Fig. H inset includes the ip(2S) region with a logarith- 
mic scale for comparison; 11802+110 '0(25') events are 
observed, consistent with the expectation of 12142+809 
ip(2S) events. We search for sources of backgrounds that 
contain a true J/0 and peak in the ir + ir~ J/ip invariant 
mass spectrum. The possibility that one or both pion 
candidates are misidentified kaons is checked by recon- 
structing the K + K~J/tp and K ± ir T J/ip final states; we 
observe featureless mass spectra. Similar studies of ISR 
events with a tt + tt~J/i/j plus one or more additional pi- 
ons reveal no structure that could feed down to produce 
a peak in the 7r + 7r~ J/0 mass spectrum. Two-photon 
events are studied directly by reversing the requirement 
on the missing mass; the number of events inferred for 
the signal region is a small fraction of those observed 
and their mass spectrum shows no structure. Hadronic 
e + e~ — * qq events pro duce J/ip at a rate that is surpris- 
ingly large [12t ll.l Il4l |l5j , but no structure is observed 
for this background. 

We evaluate the statistical significance of the enhance- 
ment using unbinned maximum likelihood fits to the 
7r+7T~ J/0 mass spectrum. To evaluate the goodness of 
fit, the fit probability is determined from the x 2 and the 
number of degrees of freedom for bin sizes of 5, 10, 20, 
40, and 50MeV/c 2 . Bins are combined with higher mass 



neighbors as needed to ensure that no bin is predicted 
to have fewer than seven entries. We try first-, second-, 
and third-order polynomials as null-hypothesis fit func- 
tions. The x 2 -probability estimates for these fits range 
from 10~ 16 to 10~ n . No substantial improvement is ob- 
tained by including 0(4040), 0(4160), or 0(4415) [ll| 
terms in the fit. We conclude that the structure near 
4.26GcV/c 2 is statistically inconsistent with a polyno- 
mial background. Henceforth, we refer to this structure 
as the F(4260). 

It is important to test the ISR-production hypothe- 
sis because the J PC = 1 assignment for the Y"(4260) 
follows from it. The ISR photon is reconstructed in 
(24+8)% of the Y~(4260) events, in agreement with the 
25% observed for ISR 0(25) events. Kinematic distri- 
butions for the signal are obtained by subtracting scaled 
distributions for events with 7r + 7r~ J/0 mass in the re- 
gions [3.86, 4.06]GeV/c 2 and [4.46, 4.66] GeV/c 2 from 
those with Tt + n~ J/0 mass in the signal region, defined as 
[4.16, 4.36] GeV/c 2 . The distribution of m 2 Rec is shown in 
Fig. [3 along with corresponding distributions for ISR 
0(25) data events and for ISR Y"(4260) Monte Carlo 
events. Good agreement is found for these distributions, 
and for all other quantities studied to test that initial- 
state radiation is responsible for these events. 

An unbinned likelihood fit to the ir + ir~ J/ip mass spec- 
trum is performed using a single relativistic Breit-Wigner 
signal function and a second-order polynomial back- 
ground. The signal function is multiplied by a phase 
space factor and convoluted with the previously de- 
scribed resolution function. The fit gives 125+23 events 
with a mass of 4259 ± 8(stat)l"g(syst) MeV/c 2 and a 
width of 88 ± 23(stat)l^(syst) MeV/c 2 . Systematic un- 
certainties include contributions from the fitting proce- 
dure, the mass scale, the mass-resolution function, and 
dependence on the model of the F(4260) — * ir + ir~ J/ip 
decay. They have been added in quadrature. Un- 
der this single-resonance hypothesis we calculate a value 
of T(y(4260) e+e") • B(Y(A260) -> n+n-J/ip) = 
5.5 ± l.Olo^eV/c 2 . The fit probability determined from 
the x 2 and the number of degrees of freedom ranges from 
0.3% to 6.6% for the same set of binning choices and 
background parameterizations used to evaluate the null 
hypothesis. To estimate the significance of the Y"(4260) 
structure conservatively, we use instead of our optimized 
selection criteria, the criteria developed in analyzing just 
the first 124fb _1 of data. Using these, we compare fits to 
the remaining 109 fb -1 of data sample with and without 
the resonance parameters determined by the first data 
sample. Using the binnings described above, we find a 
significance in the second independent data sample alone 
of 5 to 7a. The likelihood and x 2 differences between sig- 
nal and null-hypothesis fits to the full sample correspond 
to significances of at least 8a. 

The robustness of the F(4260) signal is tested with 
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FIG. 1: The -k + tt~ J/ip invariant mass spectrum in the range 
3.8—5.0 GeV/c 2 and (inset) over a wider range that includes 
the ip{2S). The points with error bars represent the selected 
data and the shaded histogram represents the scaled data 
from neighboring e + e~ and fi + \x~ mass regions (see text). 
The solid curve shows the result of the single-resonance fit 
described in the text; the dashed curve represents the back- 
ground component. 



single-resonance fits to the tt + tt~J/i/j mass spectrum for 
e+e~ and modes separately, which yield 49±16 

and 76±13 signal events, respectively. Fits give 76±18 
events for the original 124 fb _1 data set and 56±13 events 
for the next, independent 109fb _1 data set. Fitting sam- 
ples with and without reconstructed ISR photons gives 
30±11 and 96±15 events, respectively. We find consis- 
tent values for the F(4260) and the i/)(2S) when deter- 
mining the fraction of the total signal found in each of 
these subsets. 

Several additional systematic checks have been per- 
formed. Each selection criterion has been tightened (loos- 
ened) and the decrease (increase) in the signal yield is 
consistent with that for the ip(2S) data. Events selected 
when the selection criteria are reversed, individually or 
in pairs, are studied; in no case is there a significant dip 
in the signal-mass region that might indicate a bias in 
the selection procedure. 

Since the single-resonance fit probability is low we con- 
sider the possibility that the observed signal is due to two 
interfering resonances. Two-resonance fits with an inter- 
ference term find one resonance mass close to the mass 
from the single-resonance fit, but with a width as low 
as 50MeV/c 2 , plus a second narrow resonance around 
4.33GeV/c 2 . However, the fit probabilities are not sig- 
nificantly improved by two-resonance hypotheses. The 
size of our sample does not allow a statistically signifi- 
cant discrimination; we can neither exclude nor establish 
a multi-resonance hypothesis. 

The dipion invariant mass distribution for the F(4260) 
is shown in Fig. [3] Each point represents the yield of a 
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FIG. 2: The distribution of m 2 Rec . The points represent the 
data events passing all selection criteria except that on m Rec 
and having a -K + -K~J/ip mass near 4260 MeV/c 2 , minus the 
scaled distribution from neighboring -k + tt~ J/ip mass regions 
(see text) . The solid histogram represents ISR Y Monte Carlo 
events, and the dotted histogram represents the ISR tp(2S) 
data events. 
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FIG. 3: The dipion mass distribution for F(4260) -> 
tt + -k~ J/ip data is shown as points with error bars. The his- 
togram shows the distribution for Monte Carlo events where 
F(4260) — > -K + TT~J/ip is generated according to an S-wave 
phase space model. 



single-resonance fit to the tt + tt~ J/ip mass distribution for 
that 7r + 7r - mass bin. 

No enhancement has been observed in the cross section 
for e + e~ — > hadrons |ll| at energies corresponding to 
the F(4260). We compute the cross section for e + e~ — > 
■K + 7r~J/ip production at 4.25 GeV, corresponding to the 
highest bin in our data, to be about 50 pb. The inclusive 
hadronic cross section at ^/s = 4.25 GeV is 14.2 nb (llj . 
The ratio, approximately 0.34%, is smaller than the 4% 
experimental uncertainty for the hadronic cross section, 
so this mode would not have been visible. However, if 
the branching fraction of y(4260) to ir + TT~J/ip is very 
small, decays to other hadronic modes like DD would 
have been observable. This indicates that the branching 
fraction to tt + tt~ J/ip must be large compared to that for 
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^(3770) [3. 

In summary, we have used initial-state radiation events 
to study the process e+e~ — * ir + ir~ J/ip across the char- 
monium mass range. In addition to the expected ip(2S) 
events, we observe an excess of 125±23 events centered 
at a mass of ^4.26 GeV/c 2 , signifying the presence of one 
or more previously unobserved J PC = 1 states con- 
taining hidden charm. At the current level of statistics 
we are unable to distinguish the number of new states; 
the data can be characterized by a single resonance of 
mass ~4.26GeV/c 2 and of width ~90MeV/c 2 . 
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